Transition from Collisionless to Hydrodynamic Behavior in an Ultracold Atomic Gas 



S. D. Gensemer and D. S. Jin 
JILA, Quantum Physics Division, National Institute of Standards and Technology and University of Colorado, Boulder, 

Colorado 80309 



o 

o ■ 



(N 



43 

O 
CD 



■4— > 

CO 



I 

C 

O 

o 



in 
o 



13 

o 
o 



X 
S3 



Relative motion in a two-component, trapped atomic gas 
provides a sensitive probe of interactions. By studying the 
lowest frequency excitations of a two spin-state gas confined 
in a magnetic trap, we have explored the transition from the 
collisionless to the hydrodynamic regime. As a function of 
collision rate, we observe frequency shifts as large as 6% as 
well as a dramatic, non-monotonic dependence of the damping 
rate. The measurements agree qualitatively with expectations 
for behavior in the collisionless and hydrodynamic limits and 
are quantitatively compared to a classical kinetic model. 



Because dilute quantum gases have controllable inter- 
actions that can be described from first principles they 
provide an ideal testing ground for many-body theories 
of quantum fluids. Elementary excitations, in particular, 
play a key role in understanding the behavior of quan- 
tum fluids. The nature of these excitations varies dras- 
tically between two regimes, collisionless and hydrody- 
namic, based on the relative strength of interactions. In 
spite of their promise of controllable interactions, atomic 
gas experiments have had difficulty exploring this full 
range of behavior, in part because the high densities re- 
quired to reach the hydrodynamic regime also lead to 
large inelastic collision rates and consequent rapid num- 
ber loss. The recently created Bose-Einstein condensates 
in metastable He, on the other hand, are deeply in the 
hydrodynamic regime pLp| . In this work we have utilized 
relative motion in a two-component alkali atom gas to 
observe the transition from the collisionless to the hy- 
drodynamic regime. 

In general, the underlying character of excitations, as 
well as their frequencies and damping rates, depends 
on the relative strength of the interactions in the sys- 
tem. In the collisionless regime, the collision rate T co u 
is much smaller than the excitation frequency u. In 
this limit, there are few scattering events per oscillation. 
Classically, the motion is described by the single-particle 
Hamiltonian and collisions tend to damp excitations. For 
a quantum Fermi gas, this is the regime of zero sound, 
which is a collective excitation due to the self-consistent 
mean-field of a large number of particles. The oppo- 
site limit of large collision rate, T co a S> u), is called the 
hydrodynamic or collisional regime. Here, the motion 
consists of collective excitations in which the high colli- 
sion rate maintains local equilibrium throughout the gas. 
This is the regime of first sound waves, which, unlike 



single-particle excitations, are only weakly damped by 
collisions. 

For an atomic gas confined in a harmonic potential the 
excitation frequency scales as the trap frequency while 
the collision rate can also be controlled by varying the gas 
density and temperature. Excitations of a trapped clas- 
sical gas have been treated theoretically for both the col- 
lisionless and hydrodynamic regimes g-§. However, ex- 
periments have typically been in the collisionless regime 
or between the two regimes ||. In order to reach the 
hydrodynamic regime at typical densities for alkali atom 
experiments, we excite the lowest frequency excitation 
of a two-component, magnetically confined gas. The ex- 
citation mode involves center-of-mass motion along the 
weaker, axial direction of a cylindrically symmetric mag- 
netic potential. Using excitations in the weakest direc- 
tion of the trap has the advantage that the strength of 
the interactions relative to the excitation frequency can 
be increased by tightly confining the gas radially. For 
a single-component gas this dipole mode, or "slosh," is 
unaffected by interactions since collisions between atoms 
cannot alter center-of-mass momentum. However, for a 
two-component gas, consisting of atoms in two different 
spin-states for example, collisions can impact relative mo- 
tion of the two gases and this lowest frequency mode 
becomes a useful probe of interactions. Excitations of 
this type have been observed in studies of two-component 
Bose-Einstein condensates Hf[o) and have been explored 
theoretically for degenerate Fermi gases • 

For this experiment the gas consisted of 40 K atoms 
in the internal states / = 9/2, m/ = 9/2 and / = 9/2, 
rrif = 7/2 (denoted here by 1 9/2) and |7/2), respectively), 
where / is the total atomic spin and my is the magnetic 
quantum number. Because of their different magnetic 
moments, atoms in these two states have slightly differ- 
ent single-particle oscillation frequencies in the magnetic 
trap. This fact is essential to these experiments because 
we find that the crossover to the hydrodynamic regime 
for the slosh mode occurs at a collision rate set by the 
small difference in the axial trap frequencies for the two 
components. 



As described in previous work |13|, the atoms are pre- 
cooled in a vapor-cell magneto-optical trap and then 
loaded into a magnetic trap where they are further cooled 
by forced evaporation. For this experiment a gas of be- 
tween 0.35 and 3.5 million atoms, with a spin mixture 
of 45% |9/2) (55% |7/2)), was cooled to a temperature 
between 0.5 and 2 fiK. The magnetic trap strength cor- 
responded to an axial frequency of 19.84 Hz for a |9/2) 
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atom. The radial trap frequency, which was varied to ac- 
cess different relative interaction strengths in the gas, was 
set to either 135 Hz or 256 Hz. 40 K atoms are fermions, 
however the measurements described in this Letter were 
performed in the classical regime by keeping the temper- 
ature of the gas above the Fermi temperature. 

To excite a slosh mode we applied an additional mag- 
netic field that shifted the trap center along the axial 
direction. After 28 ms (roughly half of the period of os- 
cillation in the trap) this external field was switched off, 
and the motion of the atoms was allowed to evolve freely 
in the magnetic trap for time t. After this delay time, the 
trap was switched off and the gas ballistically expanded 
for 11 ms before an absorption image of the atoms was 
taken. During the ballistic expansion an inhomogeneous 
magnetic field was applied in order to spatially separate 
atoms in the two spin states through the Stern-Gcrlach 
effect [ pT|Jl^1 . Using Gaussian fits to the images recorded 
on a CCD camera we extracted the temperature, num- 
ber, and density of both components. The time evolution 
was mapped out by repeating this sequence for different 
times t and recording the position of the center-of-mass 
of each component gas after expansion. 

For a spin-polarized gas of either spin state the motion 
of the cloud center fit a sine function, with no damping 
observed over 1 sec. The axial trap frequency was thus 
determined to be cjg/27r = 19.84 Hz for the |9/2) atoms 
and ujt/2tt = 17 AA Hz for the j 7/2) atoms. In contrast, 
excitations of a two spin-state gas in the magnetic trap 
depended strongly on collisions in the gas. Because the 
atoms are fermions, s-wave collisions are forbidden be- 
tween identical atoms. Additionally, p-wave (and higher 
order) collisions are energetically forbidden at the tem- 
peratures of interest |l3[ |. However, collisions can occur 
between atoms in different spin states and it is these col- 
lisions that impact relative motion in the two-component 
system. 

Typical data are shown in Fig. |l| On the vertical axis 
are plotted the center positions of the two atom clouds, 
after 11 ms of free expansion. The data in Fig. |l](a) is 
for a small number of atoms, resulting in a relatively low 
collision rate. The |7/2) and |9/2) clouds oscillate close 
to their respective bare trap frequencies, but the motion 
is clearly damped. In Fig. 0(c), we look at the cloud 
motion in the hydrodynamic regime. Here a larger num- 
ber of atoms and tighter radial confinement is used to 
reach a much higher collision rate. The 1 7/2) and 1 9/2) 
clouds now oscillate synchronously at an intermediate fre- 
quency between the two single-particle frequencies. This 
collective mode also damps, at a rate similar to the low 
collision rate data in Fig. |l](a). Data in Fig. |l|(b) show 
the cloud motion in the transition region between the 
collisionless and hydrodynamic limits. Here the motion 
of the two clouds is coupled and the frequencies of mo- 
tion are significantly shifted from the single-particle trap 
frequencies. Furthermore, the excitation is subject to 



strong damping. 
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FIG. 1. Typical data showing the motion of the 1 9/2) (•) 
and |7/2) (o) cloud centers (z) following an abrupt shift of 
the trap potential. Traces (a) and (c) correspond to the col- 
lisionless and hydrodynamic regimes respectively, while (b) is 
in the transition region. The lines are fits to a superposition 
of two damped harmonic oscillator modes, with the solid and 
dashed lines corresponding to the motion of the |9/2) (•) and 
|7/2) (o) cloud centers, respectively. 

The lines in Fig. [l] are fits to two modes of damped, 
harmonic oscillator motion. Because the gas has two 
components the motion consists of a superposition of 
two normal modes, each with an oscillation frequency u 
and an exponential damping time r. The time-dependent 
center positions of the 1 9/2) and 1 7/2) gases, zg(t) and 
Zr (t) respectively, are simultaneously fit to the following 
function: 

z 9 (t) = A 1 e~*/ ri sin(u;if + fa) + A 2 e^/ 7 " 2 sm(uj 2 t + fc) (1) 
z 7 (t) = Bi e~ t/Tl sin(cJit + 0i) + B 2 e~ t/T2 sin(w 2 < + 2 ) 



The measured frequencies, loi/2tt and ^/^tt, for dif- 
ferent collision rates T co u are plotted in Fig. ^. The col- 
lision rate was varied primarily by changing the trapped 
gas density, either by changing the total number of atoms 
or by changing the radial trap strength. The average col- 
lision rate per atom in the trap is T co u — 2nav/N where 
N is the total number of atoms. The density overlap in- 
tegral n — J ng(r)n7(r)d 3 r and the mean relative speed 
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v were determined from two-dimensional Gaussian fits to 
the absorption images of the expanded gas. The collision 
cross-section is given by a = iwa 2 where the triplet scat- 
tering length for 40 K is a — 169a Q (a a is the Bohr radius) 

M. 



lisions are strong enough to impede single-particle exci- 
tations but are not sufficient to give rise to a well-defined 
sound wave. As expected the damping rate 1/r increases 
linearly with T in the collisionless regime (see inset to Fig. 
0) and decreases as 1/r in the hydrodynamic regime. 
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FIG. 2. Frequencies of the two excitation modes as a func- 
tion of collision rate T co u. The two modes frequencies corre- 
spond to the bare trap frequencies for the two spin-states in 
the limit of low collision rate. As F co u grows the mode fre- 
quencies shift away from the bare trap frequencies. However 
a change in behavior occurs at F co n « §(wg — L07), after which 
only one mode, with a frequency that is independent of T co u, 
is observed. Results of a classical kinetic model (lines) are 
compared to the data. 

As the collision rate was increased the measured fre- 
quencies shifted from the single-particle values. This cor- 
responds to leaving the collisionless limit where the inter- 
actions do not significantly affect the atomic motion. At 
larger r co ;/, the observed frequencies do not depend on 
collision rate. Furthermore, the frequency and damping 
time could only be extracted for one normal mode. This 
marks the transition to the regime of hydrodynamic or 
first sound where T co u is high enough to give rise to a col- 
lective excitation in the classical gas. As can be seen in 
Fig. |l](c) the observed collective mode corresponds to in- 
phase motion of nearly equal amplitude for the two com- 
ponent gases. The frequency of this mode lies between 
the two bare frequencies with a shift of approximately 
6% compared to either lo$ or oj-j. 

The measured exponential damping times, T\ and T2, 
are plotted in Fig. ^. We observe a minimum in r (or 
correspondingly a maximum in the damping rate) near 
the point where collective mode emerges. This maxi- 
mum marks the transition between the collisionless and 
hydrodynamic regimes and is well-known in the context 
of quantum fluids [0. In this transition regime the col- 



FIG. 3. Exponential damping times for the two modes of 
excitation. The damping times exhibit a minimum at the 
transition from the collisionless to the hydrodynamic regime. 
In the hydrodynamic regime, the measured r is seen to depend 
linearly on r co n . The inset shows the damping rate 1 /r which 
scales linearly with r co H in the collisionless regime. 

In Figs. ||and|^ we also compare the data to the results 
of a classical kinetic model. This model reproduces the 
previous results by Vichi and Stringari |11| but can also 
accommodate different trap frequencies for the two spin- 
states as well as different numbers of atoms in the two 
components. We model the system with the following 
coupled equations of motion: 
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u 7 z 7 , 



where F d = ±m N T coU (z g - z 7 ). 

Here Ng and N 7 are the numbers of atoms in each 
spin-state, and m is the atom mass. The viscous damp- 
ing force due to the collisional interactions in the gas, 
F4, was be derived by considering the effect of individ- 
ual collision events on the cloud's center-of-mass motion 
and integrating over all possible collisions. The expres- 
sion Fd assumes a classical gas and a small amplitude for 
the slosh. The solution of these equations consists of two 
normal modes whose frequencies and exponential damp- 
ing times are shown as the lines in Figs. |^ and |3|. For 
these theory lines the bare trap frequencies, as well as 
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the measured temperature, number, and spin composi- 
tion (45% 1 9/2)), are input to the model. In addition, we 
allow a single free parameter that is a multiplicative fac- 
tor that scales the collision rate axis. The best fit scaling 
corresponds to 1.30±0.02, which implies that the experi- 
mentally determined collision rates are low. This scaling 
factor is consistent with our estimated ±50% systematic 
uncertainty in extracting N from absorption images. 

The model gives excellent agreement with the data for 
both frequency shifts and damping times. The model re- 
veals the existence of a second, strongly damped collec- 
tive mode in the hydrodynamic regime. This collective 
mode corresponds to the spin-dipole mode studied by 
Vichi and Stringari for a degenerate Fermi gas [fill and 
is overdamped in the hydrodynamic regime. This type of 
excitation is also observed as the giant dipole resonance 
in nuclei. The model can also answer the question of 
what sets the scale for the transition from the collision- 
less to the hydrodynamic regime in the two-component 
system. Varying the bare trap frequencies in the model 
reveals that the maximum in the damping rate, which 
marks the emergence of hydrodynamic behavior, scales 
with the frequency difference of the bare modes [Q . 

Finally in Fig. |we show the emergence of coupled mo- 
tion in the two-component gas by plotting the amplitude 
ratios Ai/Bi and A2/B2 extracted from the fits to Eqn. 
1. Each normal mode is in general a linear combination 
of motion of the two component gases. A finite ratio cor- 
responds to a nontrivial combination and Fig. ^ reveals 
that the motion of the two species becomes coupled as 
T co u increases. 




40 60 

r,,, (Hz) 



FIG. 4. Ratio of the amplitudes for |9/2) and |7/2) motion. 
The ratios of the amplitudes, A\/Bi and A2/B2, are plotted 
versus collision rate V co u. As F co u increases the motion of 
the two gases becomes increasingly coupled until a collective 
mode, consisting of equal amplitude motion of the two gases, 
emerges in the hydrodynamic regime. 



In conclusion, by exploiting relative motion in a two- 
component trapped atom gas we have examined the tran- 
sition from the collisionless to the hydrodynamic regime. 
The ability to access the full range of excitation behaviors 
demonstrated here will be extremely useful for studying 
elementary excitations of quantum gases, whether Bose- 
Einstein condensates or Fermi degenerate gases. In the 
Fermi system, for example, quantum statistical reduction 
in collision rate could be revealed in a study of excita- 
tions in the hydrodynamic regime p9| . Although it has 
been previously suggested that this regime would be dif- 
ficult to reach for 40 K experiments pd| , ^9| , the results of 
this work prove otherwise. Excitations in the hydrody- 
namic regime could also be used to reveal the onset of 
superfluidity |20 2^|. In addition, spin excitations such 
as described here are interesting in their own right for 
quantum degenerate gases. 
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